home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
CD ROM Paradise Collection 4
/
CD ROM Paradise Collection 4 1995 Nov.iso
/
program
/
nrpas13.zip
/
KSONE.DEM
< prev
next >
Wrap
Text File
|
1991-04-29
|
1KB
|
54 lines
PROGRAM d13r10(input,output);
(* driver for routine KSONE *)
CONST
npts=1000;
eps=0.1;
TYPE
gldarray = ARRAY [1..npts] OF real;
glsarray = gldarray;
VAR
glinext,glinextp : integer;
glma : ARRAY [1..55] OF real;
gliset : integer;
glgset : real;
i,idum,j : integer;
d,factr,prob,varnce : real;
data : gldarray;
(*$I MODFILE.PAS *)
(*$I SORT.PAS *)
(*$I PROBKS.PAS *)
(*$I RAN3.PAS *)
(*$I GASDEV.PAS *)
(*$I ERFCC.PAS *)
FUNCTION func(x: real): real;
VAR
y : real;
BEGIN
y := x/sqrt(2.0);
func := 1.0 - erfcc(y)
END;
(*$I KSONE.PAS *)
BEGIN
gliset := 0;
idum := -5;
writeln('variance ratio':19,'k-s statistic':16,'probability':15);
writeln;
FOR i := 1 to 11 DO BEGIN
varnce := 1.0+(i-1)*eps;
factr := sqrt(varnce);
FOR j := 1 to npts DO BEGIN
data[j] := factr*abs(gasdev(idum))
END;
ksone(data,npts,d,prob);
writeln(varnce:16:6,d:16:6,prob:16:6)
END
END.